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Abstract 



We compute non-perturbatively the renormalization constants of quark bilinears on 
the lattice in the quenched approximation at three values of the coupling (3 = 6/gg = 
6.0, 6.2, 6.4 using both the Wilson and the tree-level improved SW-Clover fermion 
action. We perform a Renormalization Group analysis at the next-to-next-to-leading 
order and compute Renormalization Group invariant values for the constants. The 
results are applied to obtain a fully non-perturbative estimate of the vector and 
pseudoscalar decay constants. 
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1 Introduction 



Lattice QCD is an unique tool to compute non-perturbatively from first prin- 
ciples the mass spectrum, leptonic decay constants and in general hadronic 
matrix elements of local operators. Renormalization constants, relating the 
operators on the lattice to the continuum are necessary to extract physical 
informations from Monte Carlo simulations. In this paper we study the renor- 
malization properties of composite bilinear operators [1,2] with the quark ac- 
tion discretized a Id Wilson. 

In principle, the renormalization of quark bilinears can be computed in one- 
loop perturbation theory, as there are no power divergences [3,4]. It is well 
known, though, that lattice perturbation theory is ill behaved, due to the pres- 
ence of tadpole-like diagrams [5,6] and at values of the coupling (3 = Q/g^ = 
6.0 — 6.4, the higher-order corrections may not be small, thus introducing a 
large uncertainty in the calculation of the renormalized matrix elements in 
some continuum scheme. These problems are avoided using non-perturbative 
(NP) renormalization techniques [7,8]. The procedure proposed in [7] allows 
a full non-perturbative computation of the matrix elements of composite op- 
erators in the Regularization Independent (RI) scheme [7,9]. The matching 
between the RI scheme and MS, which is intrinsically perturbative, is com- 
puted using only continuum perturbation theory, which is well behaved. This 
method has been shown to be quite successful in reproducing results obtained 
by other methods, such as chiral Ward Identities [10]. Quite impressive is 
also the influence of the NP renormalization in the measurement of the quark 
masses [11], the chiral condensate [12] and, for four-fermion operators, in the 
restoration of the correct chiral behaviour of the .^-parameters in weak decays 
[13-15]. Moreover, any attempt to tackle the question of the AI = 1/2 rule 
must rely on NP methods [16]. 

In this paper, we extend the exploratory computations done in [7] to a high 
statistics study with both the standard Wilson action and the tree-level im- 
proved SW-Clover action [17] at different values of the coupling in the quenched 
approximation. We compare the dependence on the renormalization scale of 
the renormalized operators with respect to the solution predicted from the 
Renormalization Group Equation (RGE) at the next-to-next-to-leading-order 
(NNLO). Moreover, we use this analysis to estimate the systematics due to 
discretization errors on the renormalization constants. 

Recently, there has been much progress in the Symanzik on-shell improvement 
program [18], obtaining a non-perturbative determination of the renormaliza- 
tion constants with a fully 0(a) improved action [19,20]. An extension of this 
program, to take into account all terms of 0(a) including the ones proportional 
to the quark mass, has been proposed in [21]. Whereas the non-perturbative 
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improvement is the direction of future, it is true that Wilson and tree-level 
Clover actions are still widely used, e.g. for four-fermion operators relevant to 
weak decays (cf. for example ref. [22] and references therein). Moreover, the 
study of the renormalization properties in the chiral limit with an unimproved 
and partially improved action will provide further insight for the discretization 
effects in lattice QCD. 

The outline of this paper is as follows. In sec. 2 we review the non-perturbative 
method (NPM) proposed in ref. [7] and set the notation for the remainder of 
this work. The Renormalization Group (RG) analysis of the quark bilinears 
is outlined in sec. 3, while in sec. 4 we present the numerical results for the 
renormalization constants and discuss the systematic errors. In sec. 5 we apply 
these results to the the computation of the leptonic decay constants of the 
vector and pseudoscalar mesons. We finish with our conclusions in sec. 6. 



2 Non-perturbative renormalization method 



In this section we review the method of ref. [7], which we have used to com- 
pute non-perturbatively the renormalization constants of quark bilinears in 
the Regulatization Independent (RI) scheme [9]. The method imposes renor- 
malization conditions non-perturbatively, directly on quark and gluon Green 
functions, in a fixed gauge, with given off-shell external states of large vir- 
tually. Notice that in RI the renormalization conditions are independent of 
the regularization scheme but they depend on the external states and on the 
gauge used in the procedure. 

The renormalization scale /i 2 , determined from the virtuality of the external 
states p 2 , must satisfy the condition A QCD <C /i <C 0(1 /a), see [7]. 

We have worked in the lattice Landau gauge, defined by minimizing the func- 
tional 



Tr 



fl=l 



(1) 



The necessity to fix the gauge introduces a systematic uncertainty due to the 
existence of both continuum and lattice Gribov copies [24] and the numerical 
noise that they can generate. These effects are expected to die off at large 
virtuality and on the renormalization of two-quark operators have been found 
to be small, comparable to the statistical noise [23]. We are making the as- 
sumption that the Landau lattice gauge-fixing procedure brings gauge fixed 
lattice operators into the corresponding continuum ones as a — > [25]. 
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Let us consider a local lattice quark bilinear Or = qTq, where T is a Dirac 
matrixQ The renormalization condition is imposed on the amputated Green 
function computed between off-shell quark states of momentum p in the Lan- 
dau gauge 

A r (pa) = S q (pa)~ 1 G v {pa)S q (paY 1 (2) 



where Gr(pa) and S q (pa) are the non- amputated Green function and quark 
propagator, calculated non-perturbatively via Monte Carlo simulations [7]. 
The renormalization constant Z^ l (iia, go) of Or, in the RI scheme, is deter- 
mined by the condition 

Z^CH^MTr P r A r M| p2=At2 = 1, (3) 



and the renormalized operator is related to the bare one by Of 11 = Z^Or- In 
eq. (3) Pr is a suitable projector on the tree-level amputated Green function. 
In the case of the quark bilinears the projector is simply proportional to T*. 
Z q is the wave function renormalization which can be defined from the Ward 
Identity (WI) as [7] 



Z q (fia) = -ij^Tr 



dSjpoY 
dp 1 



(4) 



To avoid derivatives with respect to a discrete variable, we have used 

. 1 rrE^i^T^sin^a)^^) -1 



Z'(fia) 



12 



4£„=i, 4 sin (p^a) 



(5) 



which, in the Landau Gauge, differs from Z q by a finite term of order a 2 . The 
matching coefficient can be computed using continuum perturbation theory 
only, and up to order a 2 [26] 

§ = 1 - -^2 A l 2) + • • • ( 6 ) 

Z' q (47T) 2 q 1 ^ 



with, in the Landau gauge, 

Af = ^1^(3 + 22^-4^). (7) 



1 In the following, we shall denote with r = A, V, P, S the axial and vector currents 
and the pseudoscalar and scalar densities. 
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where N c is the number of colours and n/ the number of active quarks. Eqs. (3) 
and (5) define the constants Z r RI . From eq. (6) we obtain for Z^ 1 

= (l - |^Af) Z' r ^a) , (8) 



which satisfy the Ward identities at the NNLO. The matching between RI 
and MS, NDR requires continuum perturbation theory only [7,9]. Since both 
RI and MS, NDR respect chirality and the renormalized operators with the 
correct chiral behaviour are unique, we have Z^ 1 = Z^ and Zy 1 = Z v ls . For 
the same reason 

7RI 7RI 7MS 

—L= = =^L= = =2L_ (9) 

%MS ^MS £W ' v ' 



where Z m is the quark mass renormalization. In ref. [26] Z^ s /Z^ has been 
computed up to order a 2 . Using eq. (9) one can compute at the same order 
the matching coefficients of the scalar and pseudoscalar densities 



Z^) 



C (i) aKt*) 
Att (An) 2 



(10) 



where T = P, S and in the Landau gauge, 



8 (N 2 ~ 1) 

(N 2 C - 1) 
96 7V? 



(-309 + 3029 N 2 
-288 Cs - 576 N 2 ( 3 - 356 N c n f ) . 



where C3 = 1.20206 The dependence on the gauge and the external states 

of the RI scheme will cancel with the corresponding dependence of the match- 
ing coefficients in eq. (10), up to higher orders in continuum perturbative 
expansion and up to discretization errors. 



3 Renormalization Group Analysis 



The RGE expresses a general property of the Green's functions of a renormal- 
ized theory and therefore they are valid non-perturbatively. To study the RG 
properties of bilinears we work in the MS scheme and in the Landau gauge. 
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The generic, forward, renormalized two-point Green's function, computed be- 
tween quark states of virtuality p 2 obeys the RGE 



2 d 7r 
^ rf/i 2 + 2 



p 



2 9 n . . d 



_ 7r 

da, 2 



p 



0, (12) 



where the QCD /5-function and the anomalous dimension of the renormalized 
operator Or are gauge invariant to all orders in perturbation theory and are 
defined as: 



p(a s ) 2 d (a s \ ~ (a s y+ 2 

^2 (iz = -Eft z > ( 13 ) 



An dfi 2 \AnJ ~^ \An 

7r (« s ) = -2Z r -y--^r = f>« (14) 



In a continuum regularization which respects chirality the axial and vector 
currents do not get renormalized, i.e. Z A = Z v = 1, as can be easily shown 
through the Ward Identities that they satisfy. Since m(/x)P(/x), with P(p) 
the pseudoscalar density, is renormalization group invariant, the scalar and 
pseudoscalar densities have renormalization constants which obey Z$ = Zp = 
1/Z„Q Therefore one can express the anomalous dimension 7r of bilinear 
operators as a function of ■y m : 



7a = 7v = 0, (15) 

lP = 1S = -1m ■ 

To solve the RGE's in the NNLO approximation, the expansions of the (3 
function and anomalous dimension up to three loops is required. 

The running of the coupling constant a!^ s is given by 



«f\ 2 , _ 1 A In ln(g 2 ) 

An {q ' p ln(q 2 ) /3 3 InV) 

+ ^3^2) (ff h 2 Info 2 )-/?? h hi(g 2 ) + - ft) , (16) 

where q 2 = (;u/Aq^ d ) 2 . For the continuum MS scale parameter Aq<§ d at 

the NNLO, in the quenched approximation, we have used Aq^ d = 0.251 ± 
0.021 GeV [20]. The QCD /3-function is scheme independent only up to two 



2 We stress that we are studying the RGE evolution in a continuum renormalization 
scheme. 
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loops. The additional term of the expansion has been computed in the MS 
scheme in [27]: 



11 2 

f 3 c 3 ' iV r 



(17) 



ms _ 2857 3 (iV c 2 -l) 2 205 , 2 x 
1415 Ar2 ll(iV c 2 -l) 2 79 2 



54 



The mass anomalous dimension in the MS scheme up to three loops is given 
by [28,29]: 



Ira ^ : 

Im 



Im 



( 3 203 at2 5 Ar \ 
("4 + "12" c - ' 



3 203 

4 + IT 

129 129 Ar2 H413 Ar4 
N?+ Nt 



108 



/23 1177 

+n f iV c - - 12iV c C 3 - 12iV c 3 C s 

where £ is the Riemann zeta function. 



27 c J 



(18) 



The evolution of the renormalized bilinear operators is determined by eqs. (12) 
and (16). The solution can be expressed in the MS scheme in the form [30] 



(19) 



where 



(20) 



1 / q s (fj,y 

2 I 4tt 



with p t = (3,1(30 and 7? = 7 ?V (2A>). 

By using the eqs. (9) and (19), the evolution of the bilinear quark operators 
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at the NNLO in the RI scheme becomes 

o r RI (^) = ^^M, (21) 



where 

<* M M = f!g<rM. (22) 



Eqs. (21) and (22) define the evolution at NNLO in the RI scheme of the 
renormalization constants with the scale \i. In order to compare with the 
numerical NP results we define a Renormalization Group Invariant (RGI) 
constant as 



7RGI / \ _ Z'^\^a,m q a = 0) 
Z T (a)- , 



where 

= c^) (24) 

takes into account the mismatch between Z q and Z' cf. eq (6). Up to higher 
order terms in continuum perturbation theory and up to discretization errors, 
Zp- GI (a) should be independent of /1, in the region in which perturbation theory 
is valid, i.e. /i ^ 2 GeV, independent of the renormalization scheme, of the 
external states and gauge invariant. Being the continuum evolution already 
at NNLO, we assume any scale dependence to be dominated by discretization 
effects. As an estimate of this systematic error we will take the semidispersion 
of the values of the renormalization constants in the perturbative region. 



4 Non-perturbative renormalization constants 



The renormalization constants for the bilinears presented in this paper are 
obtained at three different gauge couplings g%, corresponding to (3 = = 
6.0, 6.2 and 6.4 using both the standard Wilson action and the tree-level im- 
proved SW-Clover fermion action [17]0. A summary of the parameters used 

3 We note that in our implementation of the tree-level improvement program the 
relationship Z$ = ^/Z rn is not true anymore because we have used non-local "]p- 
rotated" operators; we refer to [7,11] for details. 
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0.14143 
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Table 1 



Summary of parameters used in the non-perturbative calculation of the renormal- 
ization constants. 

in the NP calculation of the renormalization constants is presented in tab. 1. 
The errors have been obtained with the jacknife method, decimating 10 con- 
figurations at a time. The lattice scale a -1 for the different couplings has been 
determined from Mk* [34] and is shown in tab. 5. 

4-1 Axial and vector currents 

Let us consider first the axial and vector currents. Since each obeys a chiral 
Ward Identity [1] their renormalization constants are finite, i.e. scale inde- 
pendent. In figs. 1 and 2 we show Za and Zy, calculated in the RI scheme 
in the chiral limit for both the SW and Wilson actions, as a function of the 
renormalization scale //. As explained in [7], we expect to find a "window" in 
the range of values of \i in which the Z's are scale independent. Since at large 
values of \i the Z's will be more sensible to discretization errors, we expect 
the window to be wider as we approach the continuum. 

We can clearly see that these expectations are satisfied by the data shown 
in figs. 1 and 2, and also reported in tabs. 2 and 3. For the axial and vector 
currents, the RGI values given in the tables are simply the values of the con- 
stants calculated at iia ~ 1, since CjP^/z) for Za and Zy is totally negligible in 
comparison with the final error. The first error is statistical, while the second 
is the systematic error estimated as the semidispersion of the values in the 
region in which we believe perturbation theory to be reliable, i.e. /x > 2 GeV. 
In the SW case, the axial current shows a more pronounced systematic effect 
compared to the vector current, reflected in a reduced stability in fi of the 
plots. In the Wilson case, on the other hand, it's the vector current that is 
more fluctuating in the scale \i. Moreover, as far the comparison between the 
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Fig. 1. Renormalization of the axial current for (a) SW and (b) Wilson action as a 
function of fi for all couplings. The lattice spacing is determined from Mk* ■ 

SW and Wilson actions is concerned (at the same value of (3), we note that 
for the axial current the Wilson action seems more stable than the SW one, 
while for the vector current it's the SW action which shows a more pronounced 
plateau. 



In tab. 4 we also give the values of Za and Zy determined from the WI's 
[2,31,32]. Caution should be used in comparing the values obtained from the 
NPM and the WI's as the latter are sometimes available only for finite values 
of the quark mass, but since the mass dependence is expected to be weak, it 
is still a significant check. The agreement between the NPM and the WI's is 
good for the axial current in the SW case, around iia ~ 1, as already found 
in [7]. On the other hand, in the Wilson case the WI value at (3 = 6.0 is even 
larger than the NPM at f3 — 6.4, although the error quoted from the WI is 
so large 10%) that a significant comparison is not possible. For the vector 
current, we notice that in SW case the NP values are slightly higher than the 
WI although the discrepancy tends to diminish as (3 increases, whereas in the 
Wilson case, the comparison is more complicated. This is due to the well-know 
discretization errors which affect the hadronic matrix elements of the vector 
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Fig. 2. Renormalization of the vector current for (a) SW and (b) Wilson action as 
a function of \i for all couplings. The lattice spacing is determined from Mk* ■ 

current, in particular of the conserved one [31]. We expect the values of Z A and 
Zy to approach unity as move towards the continuum limit and the coupling 
goes to zero. All the data obtained with the NPM support this conclusion. 

In tab. 4 we also present the one-loop perturbative values [3,4], evaluated at 
pia — 1 in the RI scheme according to eq. (10). We present both Stardard PT 
(SPT), in which the bare coupling g% is used as the expansion parameter, and 
Boosted PT (BPT) in which the coupling used is [6] 



„ BPT 

a, = 



1 



ai att = 



J_9l 

(□) 4tt 



(25) 



where (□) is the expectation value of the plaquette. In all the perturbative 
calculations, the effects coming from the discretization of loop integrals on a 
finite lattice have been ignored. As already emphasized in [7], perturbation 
theory does not agree with the non-perturbative values, except in a number of 
limited cases, as the vector current. In the case of the pseudoscalar density, the 
failure is quite dramatic and has important phenomenological consequences on 
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the determination of the quark masses and the quark condensate, as explained 
in refs. [11]. Moreover, compared to the NPM, the perturbative values show 
a less pronounced variation with the coupling. Thus, even if the agreement 
improves as we approach the continuum, the discrepancy remains sizable even 
at the smallest coupling. 

The exploratory numerical calculations in [7] were carried out at /3 — 6.0 
and a single value of the hopping parameter k = 0.1425, corresponding to 
m q a ~ 0.07. We have checked that the results are compatible with ours within 
statistical errors. We have chosen not to separate the time and space com- 
ponents of the currents as they agree within errors. It is worth noting that 
the sizable difference between the renormalization constants of the different 
components and the large fluctuations found in [7] at large values of fj, 2 a 2 were 
not due to strong discretization effects, but to the choice of the external mo- 
menta. In fact, the components of the momenta p = (po,Pi,P2,P3) were chosen 
with large differences among the pi's, thus greatly enhancing the rotational 
symmetry breaking on the hypercubic lattice. 

Recently, the non-perturbative method of ref. [7] has also been applied to the 
renormalization of bilinears in ref. [33], in which "momentum" sources and 
sinks and translational invariance are used to reduce the statistical noise in 
the determination of the Z's. The method has been tested for Wilson fermions 
at (3 = 6.0, and the values of the Z's in the chiral limit calculated at fia ~ 1 
are Z A = 0.7807(8), Z v = 0.6822(7), Z P = 0.4357(17) and Z s = 0.6808(15). 
The errors quoted are purely statistical and are much smaller than the sys- 
tematic errors coming from the discretization effects, which are expected to be 
quite large in the Wilson case. With this in mind, we find that the agreement 
between the values of [33] and ours is resonable. 



4-2 Pseudoscalar and scalar densities 

The pseudoscalar and scalar densities differ from the axial and vector currents 
in that their renormalization is not finite but scale dependent. We show in 
figs. 3 and 4 the RGI values, computed according to eq. (23). We stress that 
the RGI values are only reliable from values of \x ^ 2 GeV, i.e. when continuum 
perturbation theory is to be trusted. The RGI value in tabs. 2 and 3 has been 
calculated from fia ~ 1. The systematic error is calculated in the same fashion 
as for the axial and vector currents. It is easily noticeable the difference in the 
behaviour of Zp and Zs as a function of [i. The former approaches a plateau 
only at relatively large values of n, while the latter presents a very clear and 
long plateau. This behaviour, also reflected in a much larger systematic error 
for the pseudoscalar, does not depend on the choice of the action. 
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1.012( 9)(12) 


0.863( 6)( 4) 


0.327( 4)(17) 


0.502( 8)( 9) 


0.652(10)(21) 



Table 2 



Non-perturbative values of Z^ l (m q = 0) with the SW action, for all couplings at 
several renormalization scales fi 2 a 2 . For the values at different scales the errors are 
statistical. The RGI values are computed from that one at fia ~ 1 according to the 
cq. (23) and the first error is statistical, the second systematic as explained in text. 
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RGI 


0.825( 6)(10) 


0.748( 5)(20) 


0.333( 2)(21) 


0.431( 4)( 4) 


0.773( 8)(41) 



Table 3 



Non-perturbative values of Z^r l (m q = 0) with the Wilson action, for all couplings 
at several renormalization scales n 2 a 2 . For the values at different scales the errors 
are statistical. The RGI values are computed from that one at [ia ~ 1 according to 
the eq. (23) and the first error is statistical, the second systematic as explained in 
text. 
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Table 4 



Perturbative values (Standard and Boosted PT) of Zp^m^ = 0) and non- 
perturbative values from the Wl with the SW and Wilson action. The perturbative 
Z's are evaluated at [x 2 a 2 = 1. 

While Zp and Z$ cannot be determined separately by imposing a Ward Iden- 
tity their ratio Zp/Zs can as it is scale independent. For this reason it can be 
treated in a similar fashion to Za and Zy and can be compared with the values 
obtained from the WL In fig. 5 we show Zp/Zs as function of fi. As expected, 
the behaviour of the ratio is dominated by the behaviour of Zp and agreement 
with the Wl determination seems to be obtained at rather high values of the 
renormalization scale, at which discretization effects should dominate. 
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Fig. 3. RGI values for the renormalization of the pseudoscalar density for (a) SW 
and (b) Wilson action as a function of fi for all couplings. The lattice spacing is 
determined from Mk*- 



5 Meson decay constants 



There are several interesting phenomenological quantities which can be ex- 
tracted from the matrix elements of quark bilinears, such as leptonic decay 
constants and quark masses. The problem of a non-perturbative measurement 
of quark masses has been addressed in [11] to which we refer the reader for all 
details. A very interesting by-product of both quark masses and decay con- 
stants is the estimate of the chiral condensate, which is of great phenomeno- 
logical relevance. The issue will be addressed in a forthcoming paper [12], 
together with a detailed theoretical analysis. 



In the remainder of this section, we concentrate on the determination of the 
leptonic decay constants of mesons. The pseudoscalar and vector decay con- 
stants, fps and fv, are defined as 
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Fig. 4. RGI values for the renormalization of the scalar density for (a) SW and (b) 
Wilson action as a function of /x for all couplings. The lattice spacing is determined 
from Mk*- 



(0\A \PS)=i^M PS , 

A 



(0\V t \V,r) 



Ml 



fvZ 



(26) 
(27) 



v 



where e\ is the vector-meson polarization, M PS and My are the pseudoscalar 
and vector masses, A and Vi the temporal and spatial components of the axial 
and vector currents respectively, and Z v ,a the corresponding renormalization 
constants. 



In tab. 5 we summarize the parameters used in the simulations. The bare 
lattice decay constants have been extracted from the appropriate correlation 
functions as described in ref. [34], to which we refer for details. In tab. 6 
and 7 we present the results for the decay constants, both the bare unrenor- 
malized values in lattice units and the renormalized ones in physical units, 
obtained using the renormalization constants at fia ~ 1. In the estimate of 
the error on renormalized decay constants we have neglected the errors on the 
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Fig. 5. Renormalization of the pseudo-scalar to scalar density ratio for (a) SW and (b) 
Wilson action as a function of /x for all couplings. The lattice spacing is determined 
from Mk*- 

renormalization constants. The non-perturbatively renormalized values and 
the experimental values present roughly a 10 — 15% discrepancy, which is 
to be expected, considering that we have performed all calculations in the 
quenched approximation. It is important, though, to note that compared to 
the values obtained with perturbative renormalization in ref. [34], the values 
obtained with a non-perturbative renormalization are in general closer to the 
experimental results. We have not attempted an extrapolation to the con- 
tinuum limit as the physical volume at the smallest coupling is too small to 
confidently extract hadronic matrix elements. 



6 Conclusions 



In this paper we have performed a systematic study of the renormalization of 
quark bilinears, in a non-perturbative fashion. We have also analyzed the dis- 
cretization effects, by performing our calculations with two different actions, 
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Table 5 

Summary of the parameters used in the calculation of the matrix elements. 



Action 


6.0 6.0 6.2 6.2 6.4 6.4 
SW Wilson SW Wilson SW Wilson 


(fx a)/Z A 
(/tt a)/Z A 


0.0735(18) 0.0944(26) 0.0540(23) 0.0640(21) 0.0406(14) 0.0480(16) 
0.0661(21) 0.0878(31) 0.0470(28) 0.0568(24) 0.0370(16) 0.0438(19) 


1/Uk*Z v ) 
V(fpZv) 


0.348(9) 0.451(11) 0.332(8) 0.417(8) 0.284(9) 0.363(8) 
0.366(14) 0.482(16) 0.359(20) 0.446(15) 0.290(13) 0.377(12) 
0.384(20) 0.513(21) 0.386(33) 0.475(22) 0.297(17) 0.391(17) 



Table 6 

Lattice bare decay constants for all couplings and both actions. 



P Exp 
Action 


6.0 6.0 6.2 6.2 6.4 6.4 
SW Wilson SW Wilson SW Wilson 


f K 0.1598 
U 0.1307 


0.163(12) 0.172(6) 0.150(8) 0.155(5) 0.164(9) 0.164(8) 
0.147(12) 0.160(7) 0.131(10) 0.138(6) 0.150(10) 0.150(10) 


i/U °- 23 

VfK* 

l/U 0.28 


0.294(8) 0.319(8) 0.281(7) 0.302(6) 0.245(8) 0.271(6) 
0.309(12) 0.341(11) 0.304(17) 0.323(11) 0.251(11) 0.282(9) 
0.324(17) 0.363(15) 0.327(28) 0.344(16) 0.256(14) 0.293(12) 



Table 7 



Non-perturbatively renormalized decay constants (in GeV) for all couplings and 
both actions. 

the standard Wilson action and the tree-level improved SW-Clover action, at 
three different values of the couplings. We have performed a RG analysis at 
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the NNLO and denned RGI values for the scale dependent renormalization 
constants. We have also used this approach to estimate the systematic error 
induced by discretization. Finally, we have applied our results to the calcu- 
lation of the pseudoscalar and vector decay constants and we find that the 
non-perturbatively renormalized values, albeit with still sizable statistical er- 
rors, show a trend towards the experimental values with respect to the ones 
obtained with perturbative renormalization. With our data an extrapolation 
to the continuum limit is not reliable as the physical volume at the smallest 
coupling is too small to confidently extract hadronic matrix elements. 



Acknowledgements 

We warmly thank V. Lubicz, G. Martinelli, M. Testa and A. Vladikas for 
enlightening discussions. MT acknowledges the support of PPARC through 
grant GR/L22744. 



19 



References 



[I] M. Bochicchio, L. Maiani, G. Martinelli, G.C. Rossi and M. Testa, Nucl. 
Phys. B 262 (1985) 331. 

[2] L. Maiani and G. Martinelli, Phys. Lett. B 178 (1986) 265. 

[3] G. Martinelli and Y.C. Zhang, Phys. Lett. B 123 (1983) 433. 
B. Meyer and C. Smith, Phys. Lett. B 123 (1983) 62. 

[4] A. Borrelli, C. Pittori, R. Frezzotti and E. Gabrielli, Nucl. Phys. B 409 
(1993) 382. 

[5] G. Parisi, in "High Energy Physics - 1980", Proceedings of the 
XXth International Conference, Madison, Wisconsin, eds. L. Durand and 
L. G. Pondrom (American Institute of Physics, New York, 1981). 

[6] G.P. Lepage and PB. Mackenzie, Phys. Rev. D48 (1993) 2250. 

[7] G. Martinelli, C. Pittori, C.T. Sachrajda, M. Testa and A. Vladikas, Nucl. 
Phys. B 445 (1995) 81. 

[8] K. Jansen et al, Phys. Lett. B 372 (1996) 275. 

[9] M. Ciuchini, E. Franco, G. Martinelli, L. Reina, L. Silvestrini, Z. 
Phys. C 68 (1995) 239. 

[10] G. Martinelli, S. Petrarca, C.T. Sachrajda and A. Vladikas, Phys. 
Lett. B 311 (1993) 241, E: B317 (1993) 660. 

[II] C. R. Anton et al, Nucl. Phys. B431 (1994) 667. 

V. Gimenez, L. Giusti, F. Rapuano and M. Talevi, hep- lat / 9 8 1 2 1 

[12] L. Giusti, F. Rapuano , M. Talevi, A. Vladikas, in preparation. 

[13] A. Donini, G. Martinelli, C.T. Sachrajda, M. Talevi and A. Vladikas, Phys. 
Lett. B360 (1995) 83. 

[14] L. Conti, A. Donini, V. Gimenez, G. Martinelli, M. Talevi and A. Vladikas, 
Phys. Lett. B421 (1998) 273. 

[15] S. Aoki et al. (JLQCD Collaboration), Nucl. Phys. (Proc. Suppl.) 
B 60A (1998) 67. 

[16] C. Dawson, G. Martinelli, G.C. Rossi, C.T. Sachrajda, S. Sharpe, M. Talevi 
and M. Testa, Nucl. Phys. B 514 (1998) 313. 

[17] B. Sheikholeslami and R. Wohlert, Nucl. Phys. B 259 (1985) 572. 

G. Heatlie, G. Martinelli, C. Pittori, G.C. Rossi and C.T. Sachrajda, Nucl. 
Phys. B 352 (1991) 266. 



20 



[18] K. Symanzik, Nucl. Phys. B 226 (1983) 187 and 205. 

M. Liischer and P. Weisz, Commun. Math. Phys. 97 (1985) 59, E: Commun. 
Math. Phys. 98 (1985) 433. 

[19] M. Liischer, S. Sint, R. Sommer and P. Weisz, Nucl. Phys. B 478 (1996) 365. 

M. Liischer, S. Sint, R. Sommer, P. Weisz and U. Wolff, Nucl. Phys. B 491 
(1997) 323. 

M. Liischer, S. Sint, R. Sommer and H. Wittig, Nucl. Phys. B 491 (1997) 344. 

[20] S. Capitani, M. Guagnelli, M. Liischer, S. Sint, R. Sommer, P. Weisz and 
H. Wittig, Nucl. Phys. B (Proc. Suppl.) 63 (1998) 153. 

[21] G. Martinelli, G.C. Rossi, C.T. Sachrajda, S. Sharpe, M. Talevi and M. Testa. 
Phys. Lett. B 411 (1997) 141. 

[22] S. Sharpe, Nucl. Phys. B (Proc. Suppl.) 53 (1997) 181. 

[23] M. L. Paciello, S. Petrarca, B. Taglienti and A. Vladikas, Phys. Lett. B 341 
(1994) 187. 

[24] C. Parrinello, S. Petrarca and A. Vladikas, Phys. Lett. B 268 (1991) 236. 

[25] L. Giusti, Nucl. Phys. B 498 (1997) 331. 

L. Giusti, M.L. Paciello, S. Petrarca, B. Taglienti and M. Testa, 
|hep-lat/9803021 



[26] E. Franco and V. Lubicz, |hep-ph/980349l 



[27] T. van Ritbergen, J.A.M. Vermaseren and S.A. Larin, Phys. Lett. B400 (1997) 
379. 

[28] J.A.M. Vermaseren, S.A. Larin and T. van Ritbergen, Phys. Lett. B405 (1997) 
327. 

[29] K.G. Chetyrkin, Phys. Lett. B404 (1997) 161. 
[30] R. Tarrach, Nucl. Phys. B183 (1981) 384. 



[31] M. Crisafulli, V. Lubicz and A. Vladikas, hep-lat / 9707025 , to appear in 
Eur. Phys. J. C. 

[32] D.S. Henty et al. (UKQCD Collaboration), Phys. Rev. D51 (1995) 5323. 

[33] M. Gockeler et al, Nucl. Phys. B (Proc. Suppl.) 63 (1998) 868. 
H. Oelrich, private communication. 

[34] C. R. Allton, V. Gimenez, L. Giusti and F. Rapuano, Nucl. 
Phys. B489 (1997) 427. 



21 



